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Abstract 

In a recent paper (arXiv:cond-mat/0911.2514), one of us (FYW) considered the Potts model 
and bond and site percolation on two general classes of two-dimensional lattices, the triangular- 
type and kagome-type lattices, and obtained closed-form expressions for the critical frontier with 
applications to various lattice models. For the triangular-type lattices Wu's result is exact, and 
for the kagome-type lattices Wu's expression is under a homogeneity assumption. The purpose 
of the present paper is two-fold: First, an essential step in Wu's analysis is the derivation of 
lattice-dependent constants A, B, C for various lattice models, a process which can be tedious. 
We present here a derivation of these constants for subnet networks using a computer algorithm. 
Secondly, by means of a finite-size scaling analysis based on numerical transfer matrix calculations, 
we deduce critical properties and critical thresholds of various models and assess the accuracy of the 
homogeneity assumption. Specifically, we analyze the (/-state Potts model and the bond percolation 
on the 3-12 and kagome-type subnet lattices (n x n) : (n x n), n < 4, for which the exact solution 
is not known. To calibrate the accuracy of the finite-size procedure, we apply the same numerical 
analysis to models for which the exact critical frontiers are known. The comparison of numerical 
and exact results shows that our numerical determination of critical thresholds is accurate to 7 or 8 
significant digits. This in turn infers that the homogeneity assumption determines critical frontiers 
with an accuracy of 5 decimal places or higher. Finally, we also obtained the exact percolation 
thresholds for site percolation on kagome-type subnet lattices (1 x 1) : (n x n) for 1 < n < 6. 

PACS numbers: 05.50.+q, 64.60. Cn, 64.60.Fr, 75.10.Hk 
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I. INTRODUCTION 

The g-state Potts model jl, 2] is a very important model in the study of phase transitions 
and critical phenomena. The critical frontier, or the loci of critical points, of the Potts model 
was first determined by Potts [l| for the square lattice. The critical exponents of the Potts 
model are obtained by conjectures on the basis of numerical evidence and by using Coulomb 
gas theory js-7|. According to the universality hypothesis jsj, the Potts model on different 
lattices belongs to the same universality class. But the determination of the critical frontier 
of the g-state Potts model in general, which includes the q — 1 bond and site percolation, is 
still an outstanding challenge. Particularly, the threshold of site percolation has remained 
largely unresolved. 



In a recent paper |9|, hereafter referred to as I, one of us (FYW) considered the Potts 
model on two classes of very general two-dimensional lattices, the triangular-type and 
kagome-type lattices shown in Fig. [TJ The Boltzmann weights W of the hatched trian- 
gles denote interactions involving 3 spins ri, r 2 , t 3 = 1,2, ... ,q surrounding a triangle, and 
are given by 

Wa(1,2,3) = A + B(5 12 + 5 23 + 5 31 ) + C5 123 , 

W v (l, 2, 3) = A' + B'(S 12 + 5 23 + 5 31 ) + C'5 123 , (1) 

where Sij = #Kr(si, s j)> ^123 = ^12^23^31, and A, B,C, A', B r ,C are constants. Spin interac- 
tions within the hatched areas can be either 2- or 3-site couplings. The hatched triangles 
can have internal structures such as the stack-of-triangle subnets, which are of recent in- 
terest 
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11], shown in Fig. [21 We refer to these structures as subnet networks. These 
stack-of-triangle lattices are called subnet lattices. Examples of triangular subnet lattices 
and kagome-type subnet lattices are shown in Fig. 2 and 3 of I. The lxl subnet lattices 
are the triangular and kagome lattices themselves. We shall call a kagome-type lattice with 
m x m down-pointing and n x n up-pointing subnets an (m x m) : (n x n) subnet lattice. 

In I, Wu derived closed-form expressions for the critical frontier of the g-state Potts model 
for the 2 types of lattices in Fig. [TJ For the triangular-type lattices the critical frontier is 
exact, but for the kagome-type lattices the critical frontier is obtained under a homogeneity 
assumption. 

The purpose of this paper is two-fold: 
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(a) 



(b) 



FIG. 1: (a) Triangular-type lattice, (b) Kagome-type lattice. 




(a) (b) 



(c) 



(d) 



FIG. 2: Subnet networks, (a) 1 x 1 subnet, (b) 2 x 2 subnet, (c) 3 x 3 subnet, (d) 4 x 4 subnet. 
The dots denote triangles with 3-site interactions. 

First, an essential step in Wu's analysis is the derivation of relevant lattice-dependent 
constants A,B,C for subnet networks. The derivation, while elementary, is tedious. Here 
we use a computer algorithm to evaluate them. Details of the algorithm are described in 



Secondly, we determine the critical frontier numerically and examine the accuracy of 
the homogeneity assumption. Specifically, we carry out a finite-size scaling analysis based 
on transfer matrix calculations to numerically determine the critical frontier for several 
lattice models, including the Potts model on the 3-12 and kagome-type (n x n) : (n x n) 
subnet lattices, for which the exact thresholds are not known. To assess the accuracy of the 
numerical determination, we also apply the procedure to models for which the exact critical 
thresholds are known. These include the Ising model and site percolation on the 3-12 lattice 
and kagome-type (1 x 1) : (n x n) subnet lattices, n < 6. Comparison of numerical and 
known exact results shows that our numerical procedure is accurate to 7 or 8 decimal places. 
This in turn infers that the critical frontier determined using the homogeneity assumption 
(qI of I is accurate to 5 decimal places or higher. 



Sec. EH 
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Our paper is organized as follows: The main findings of I are summarized in Sec. [TT1 
We describe in Sec. II III the algorithm we use to obtain the expressions of A, B, C for the 
Potts model with pure 2- and/or 3-site subnet interactions. The resulting expressions of 
A, B, C are listed in the Appendix. In Sec. IIV[ we describe the transfer matrix technique 
and the finite-size scaling method. Numerical results of our transfer matrix calculations and 
finite-size scaling analysis are given in Sec. |VJ New exact thresholds are also given in Sec. 
IVl for site percolation on kagome-type (1 x 1) : (n x n) subnet lattices for n up to 6. We 
summarize our main findings in Sec. I VI I 



II. MAIN RESULTS OF I 



We summarize in this section the main results of I. 

For the triangular-type lattice shown in Fig. [IJa), the partition function is 

Z^; A,B,C) = J2U WaV, J, k) , (2) 
n=i A 



where the products are taken over the up-pointing triangles. Wu [s[ Q showed that, in the 
regime 

2B + C>0, 3B + C>0, (3) 
in which the ground state of is ferromagnetic, the exact critical frontier is given by 

qA = C. (4) 

The critical function (j3J) yields the exact thresholds of site percolation on lattices gen- 
erated from triangular-type lattices. Consider a Potts model on a triangular subnet lattice 
with pure 3-site interactions in doted triangles shown in Fig. |2J Regarding faces of 3-spin 
interactions in an n x n subnet as sites of a new lattice, the Potts model maps to a site 
percolation on a (1 x 1) : (n — 1) x (n — 1) kagome-type subnet lattice. The critical frontier 
(jlj) then gives the exact threshold of the site percolation. Examples of the mapping are 
shown in Fig. [3] and H] for n = 2 and 3. 

For kagome-type lattices shown in Fig. QJb) the partition function is 

Z^ g (q;A,B,C;A',B',C) = [U W ^^ k ^ ' HI W v (i',f,k')] . (5) 

Ti=l A v 




FIG. 3: Site percolation on the kagome lattice. 




FIG. 4: Site percolation on the (1 x 1) : (2 x 2) kagome-type subnet lattice. 

Wu [9] obtained its critical frontier 

{q 2 A + 3qB + C)(q 2 A' + 3qB' + C) - 3{qB + C)(qB' + C) - (q - 2)CC' = , (6) 

under a homogeneity assumption. 

The critical point in the case of q = 2 computed from ([6]) is exact. Wu [9] also used ([6]) 
to compute Potts thresholds for the 3-12 and the (m x m) : (n x n) kagome-type subnet 
lattices for m, n < 4 for which the exact thresholds are not known. In addition, Wu deduced 
the known exact threshold of site percolation on the 3-12 lattice by considering the Potts 
model on the (2 x 2) : (2 x 2) kagome-type lattice as shown in Fig. [51 In this case the 
homogeneity assumption turns out to give the exactly known critical frontier. 




FIG. 5: Site percolation on the 3-12 lattice. 
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III. EVALUATION OF A, B, C FOR SUBNET NETWORKS 



In this section we describe the computer algorithm we use to evaluate expressions of 
A,B,C for the Potts model with 2- and/or 3-site interactions in subnet networks. 

For the Potts model with 2-site interactions, consider the 2x2 subnet network in Fig. 
[6]^a) as an example. The Boltzmann weight is 

w A = + ^X 1 + ^X 1 + vS 2A) • • • , (7) 

{4,5,6} 

where v — e K — 1, K is the 2-site coupling of the Potts model. 

Terms in the expansion of the products can be represented by graphs. As shown in Fig. 
El^a), there are 9 bonds in the subnet. Define two states for each bond, occupied and vacant, 
then there are total 2 9 = 512 graphs corresponding to the 512 terms in the expansion of (JTj). 
For example, Fig. [6](b) is a graph that corresponds to the term ^1,6^2,6 = Q 2 v 2 Si j2 

{4,5,6} 

contributing to B with a term q 2 v 2 . 

The 512 graphs are divided into five types according to following rules: 

1. Type-1, graphs with isolated spins 1, 2 and 3. The sum of these graphs generates the 
expression of A. 

2. Type-2, graphs with spins 1 and 2 connected and spin 3 isolated. The sum of these 
graphs contributes to the expression of B5u- For clarity we denote it as Bw 

3. Type-3, graphs with spins 2 and 3 connected and spin 1 isolated. The sum of these 
graphs contributes to the expression of -8^23 and denoted as -823- 

4. Type-4, graphs with spins 3 and 1 connected and spin 2 isolated. The sum of these 
graphs contributes to the expression of B5 3 i and denoted as B 31 . 




FIG. 6: (a) 2 x 2 subnet with pure 2-site interactions, (b) Red (bold) bonds are occupied, other 
bonds are vacant. 
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5. Type-5, graphs with all three spins 1, 2 and 3 connected. The sum of these graphs 
gives rise to the expression of C. 

For the Potts model with uniform and symmetric interactions, we have Bi 2 = B 2 s = 
B 31 = B. 

The algorithm of our program is to generate the 512 graphs one by one, compute the 
weight of each graph, and classify them into the five types. The graph weight assumes the 
form g nc ti n « ; where n v is the number of occupied bonds in the graph, and n c is the number 
of independent clusters isolated from, i.e., not connected to, sites 1, 2 or 3. For example, 
the graph in Fig. El^b) has n v = n c = 2 and the weight q 2 v 2 . 

The algorithm of our program is therefore as follows: 

1. Generate one term, i.e., a graph, by choosing a set of occupied bonds. 

2. Count the number of independent clusters isolated from site 1, 2 or 3 as n c . 

3. Count the number of occupied bonds n v . 

4. Assign a term q nc v nv to A, B, or C according to the aforementioned rules. 

5. Go to 1 for another graph until all 512 graphs are exhausted. 

The procedure for the Potts model with pure 3-site interactions M is similar. Take the 
case shown in Fig. [7(a) as an example. The doted up-pointing triangles have pure 3-site 
interactions and the Boltzmann weight of the 2x2 subnet can be written as 




FIG. 7: (a) 2 x 2 subnet with pure 3-site interactions indicated by dots, (b) Solid dot denotes the 
dot is occupied; open dot denotes it is vacant. 
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where m = e — 1. 

To obtain an expansion in the form of (jSJ), we define two states of the dots as either 
occupied or vacant. Thus there are 2 3 = 8 graphs corresponding to the 8 terms in flS]). 
However, up to this point, clusters are defined by the connectivity of Potts sites, not by the 
dotted faces. But the connectivity can be readily translated to that of the dotted faces. A 
moment's reflections shows that the weight contributing to A, B or C is simply q nc m nm , 
where n c is the number of independent clusters not containing sites 1, 2, or 3, and n m is the 
number of occupied dots. 

The rules to divide the graphs into five types corresponding to A,B,C are the same as 
the ones for pure 2-site interactions. For example, the graph in Fig. E^b) has n m = 2, n c = 

and corresponds to the term ^2 (mSi^) (77182^) = m 2 Si 2 , thus contributing to B 12 with a 

4,5,6 

term m 2 . 

The algorithm to obtain expressions of A,B,C is therefore very similar to the one de- 
scribed in the above for 2-site interactions: 

1. Generate one term, i.e., a graph, by choosing a set of occupied dots. 

2. Count the number of clusters isolated from sites 1, 2 or 3 as n c . 

3. Count the number of occupied dots n m . 

4. Assign q nc m nm to A, B or C respectively according to the aforementioned rules. 

5. Go to 1 for another graph until all possible graphs are exhausted. 

In the Appendix we present expressions of A, B, C for the Potts model on n x n subnets 
with 2-site interactions for n < 4, and for subnets with 3-site interactions for n < 7. 

IV. THE TRANSFER MATRIX AND FINITE-SIZE SCALING 

We use the method of transfer matrix to calculate statistical variables for lattice models 
wrapped on a cylinder with circumference L and length N. For lattices shown in Fig. [1] 
with hatched triangles, L and N count up- and down-pointing hatched triangles (rather 
than individual Potts spins within each triangle). Thus, for an (m x m) : (n x n) lattice 
shown in Fig. E^a), there are actually (m + n)L Potts spins in a length L. 
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For the Potts model, we build the transfer matrix by using the random-cluster represen- 



tation of the Potts partition function 
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Z = J2v nbi9) q nc{9 \ (9) 
a 

where the summation is over all subgraphs g of the lattice (or graph) on which the Potts 
model is defined, v — e K — 1, and rib(g) and n c (g) are, respectively, the number of bonds 
and clusters in g. For the (q = 1) bond percolation we have simply Z = (1 + v ) E , where E 
is the total number of edges of the lattice. 

The concept of connectivity plays an essential role in the building of the transfer matrix. 
Sites that belong to the same cluster are said to be connected. In the N x L cylinder, each 
of the L end sites of the cylinder is either isolated from or connected to other end sites. The 
connectivity of the L end sites of the cylinder is described by non-crossing partitions of the 
sites. There are a total of 

^^( 2L r) do) 



2L + 1 V L 



such non-crossing partitions 15| indexed by (3. In the transfer matrix consideration the non- 



crossing partitions are mapped onto and coded by a set of integers 1, 2, • • • , di- A detailed 



explanation of the coding procedure can be found in 16] 



The partition function of the Potts model can therefore be written as 

zw=5>r. (ii) 

where is the partition sum restricted to the partition (3. The restricted sums and 
Zp N ~ 1S> are connected by a transfer matrix T in the form of 

^=Evr' ( i2 ) 

where 

T a ,p = J2 vAn " {9) q Anc{9) (13) 

9 

are elements of T. Clearly, T has the dimension d^x d^. It is also clear that the summation 
in (TIB"]) is over subgraphs g connecting partitions (3 and a of the (N — l)-th and iV-th rows, 
with Arib(g) and An c (g) denoting, respectively, the net (positive or negative) change of the 
number of bonds and clusters due to the introduction of g. 



To conserve computer memory and running time, the transfer matrix is converted into a 
product of sparse matrices as described below (see [16| for further details). This technique 
has proved to be very efficient in the transfer matrix study of the Potts model, the O(n) 



loop and other lattice models 



17rl23| 



The transfer matrix can be regarded as adding a new layer to the system. This process 
converts the transfer matrix T into a product of L+l sparse matrices for kagome-type lattices 
of Fig. []Jb). The first sparse matrix T\ adds one down-triangle with two 'new' corner sites 
on top of an 'old' site followed by a shift of labeling of sites. This gives rise to a new layer 
with L + l sites as shown in Figs. E(a) and[E(b). The matrix T\ is a di+i x di rectangular 
matrix. The second sparse matrix T 2 adds one up- and one down-triangle simultaneously. 
By shifting the labels of the sites of the top (new) layer cyclically, T 2 brings in two new 
corner sites L and L + 1 on top of the current layer, and covers two old corner sites L + l 
and 1 , as shown in Figs. Mb) and[8](c). T 2 is a d^ + i x <i L+1 square matrix. After L — 1 such 
steps, the graph is transformed to the one shown in Fig. [8(d). By adding an up-triangle on 
the two old sites L + l and 1 under cylindrical boundary condition, the last sparse matrix 
T3 adds the last new corner site L to the system as shown from Figs. El(d) and[H(e). Labels 
of the top sites are shifted and T 3 is a d^ x d^+i rectangular matrix. 

The transfer matrix T now assumes the form of a product of sparse matrices, 

T = T{T^- X T X . (14) 

In the actual calculation, we need to store only the positions and values of the nonzero 
elements of each sparse matrix in a few one-dimensional arrays. 

In constructing these sparse matrices, one needs to enumerate all possible graphs inside 
the added (one or two) hatched triangles for a given connectedness of the partitions of the 
new and old layers. For example, in the construction of T\, we need to add a subnet shown 
in Fig. [2] (flipped vertically). In the case of 1 x 1 subnet, it is straightforward to enumerate 
all graphs manually. However, in the case of 2 x 2 and higher order subnets, it is tedious 
and sometimes impossible to count all possible graphs by hand. Therefore, we make use of 
a computer algorithm similar to the one used in obtaining expressions of A, B, C to count 
Ant, and An c . 
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L L+l 

(a) (b) 



L L 2 3 




(e) 

FIG. 8: Construction of the transfer matrix for kagome-type subnet lattices. The procedure consists 
of several steps each corresponding to a sparse transfer matrix. 

For site percolation, the partition function is 

9 

= 1, (15) 

where the summation is over all site percolation configurations g, s is the probability that 
a site is occupied, n s (g) is the number of occupied sites in g, and N is the total number 
of sites. The corresponding transfer matrix is defined in a way similar to that of the Potts 
model, but with a twist due to the presence of vacant sites and henceforth not all L end 
sites are occupied. Denote the number of end sites that are occupied by n — 0, 1, • • • , L 
which can be distributed in ( L ) different ways. Then there are a total of 

di=j:(% n He) 

n=0 ^ ' 

non-crossing partitions and the transfer matrix has the dimension df x df. It is clear that 
we have df > d L . 
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The partitions can again be coded by means of a sequence of integers 1, 2, • • • , df . The 
coding al gori thm is the same as that used in the consideration of the Potts model with 
vacancies 18j] and in the study of site percolation 24 1. 

To determine the critical threshold of the Potts model and/or site percolation, we calcu- 
late the magnetic scaled gap 

X h (v,L) = —±- (17) 

where £h(v,L) is the magnetic correlation length (with v replaced by s for site percolation). 
In the language of the random cluster model and site percolation, the magnetic correlation 
function is defined to be the probability that two sites at a distance r belong to the same 
cluster, or 

9r = Y > 

where Z' = J2' g v nb ^q nc ^ for the random cluster model and Z' = J2' g s ns(9) (l — s) N ~ n "^ 
for site percolation. The summations in Z' are the same as in fl9]) and ( TT5|) but restricted to 
subgraphs g with at least one cluster spanning from row 1 to row r. 

We define a transfer matrix, hereafter referred to as the magnetic sector of the transfer 
matrix, based on Z', in a way similar to that of the transfer matrix based on Z in ( TT5j) 
or (115p in the 'non-magnetic' sector. In constructing the magnetic sector of the transfer 
matrix, we use the 'magnetic' type connectivity of the L end sites of the cylinder, which, 
in addition to describing how sites are connected, specifies which sites are still connected to 
a site in row 1. These sites are called 'magnetic sites'. To count the total number of non- 
crossing partitions describing the 'magnetic' type connectivity, we first code the positions of 
the magnetic sites by means of a binary number m = 0, 1, 2 L — 1, where the binary digit 1 
denotes a magnetic site. The magnetic sites divide the remaining sites in g(m) groups such 
that two sites in different groups cannot be connected. Let n(j) be the number of sites in 
the j-th group. Then there are 

h m = d n (i)d n (2) ■ ■ ■ rf n (g(m)) (19) 

non-crossing partitions for the Potts model and 

= dn(l)dn(2) ' ' ' dn(g(m)) (20) 

non-crossing partitions for site percolation. The total number of non-crossing partitions is 
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therefore 

jM _ 

m=0 



4 m) = E ^ ( 21 ) 

for the Potts model, and 

d? m, = Z>£ ( 22 ) 

m=0 

for the site percolation. 

The partitions can again be coded by means of a sequence of integers. A detailed de- 
scription of the coding algorithm can be found in [181 ] . The magnetic sectors of the transfer 
matrix now have the dimensions d™ x d^f 1 and d^} m ^ x d^} m ^ for the Potts model and 
site percolation, respectively, and are much larger than those of the non-magnetic sectors. 
The magnetic sector of the transfer matrix can also be converted into a product of sparse 
matrices in the same way as in the case of the non-magnetic sector. 

The inverse magnetic correlation length is given by 

£ h (v,L) \X' J 

where Ao and X' are the leading eigenvalues of the transfer matrix in the non-magnetic and 
magnetic sector respectively, ( is a geometrical factor which is the ratio between the unit of 
L and the thickness of a layer added by the transfer matrix. The magnetic scaled gap then 
follows. 

According to finite-size scaling theory 25] and Cardy's conformal mapping 26], Xh(v, L) 
can be expanded as 

X h (v, L) = X h + atL yt + buL y " + ... , (24) 

where is the magnetic scaling dimension, t is the deviation from the critical point, and 
u the irrelevant field. Here, y t is the thermal renormalization exponent, y u the leading 
irrelevant renormalization exponent, and a and b are unknown constants. 

We substitute ff24l) into the finite-size scaling equation connecting lattices of sizes L and 
L-l, 

X h (v,L)=X h (v,L-l), (25) 
and denote the solution of (1251) by v c (L), which has the expansion 

Vc (L) =v c + a'uL yu - yt + ■■■ , (26) 
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where a' is an unknown constant. Because y u < and y t > 0, v c (L) for a sequence of 
increasing system sizes converge to the critical point v c . 

At v c (L), the expression X h (v c (L), L) for a sequence of sizes L converge to the magnetic 
scaling dimension as 

X h (v c (L), L)=X h + b'uU>- + ... , (27) 

with b' an unknown constant. This determines the magnetic scaling dimension X^. 
The free energy per unit distance is given by 

/(i) = ^, (28) 
where Aq is the largest eigenvalue of the transfer matrix in the non-magnetic sector. Accord- 



ing to conformal invariance theory, the 



energy density at the critical point is 



arge-L asymptotic finite-size behavior of the free 



/(L) ~ /(oo) + ^ , (29) 

where c is the conformal anomaly. 

The conformal anomaly c and the magnetic scaling dimension Xh are two important 
universal quantities defining the universality class. For the two-dimensional (/-state Potts 
model, they are given by the conformal invariance theory and Coulomb gas method (3I.I4I.I27- 



29j as 



6(1 -g) 2 g 3 

°=l--44-. X k =l-»---, (30) 

where 

v ^=-2cos(7T( 7 ), \<g<l- (31) 



V. NUMERICAL AND SOME EXACT RESULTS 

In this section we present numerical results of our transfer matrix calculations and finite- 
size scaling analysis for the 3-12 and kagome-type lattices. We also present some exact 
results for site percolation on the (1 x 1) : (n x n) lattices. 

A. The g-state Potts model on the 3-12 and kagome-type subnet lattices 

Critical points are estimated by extrapolating the solutions of (|25|) for a sequence of in- 
creasing system sizes in accordance with the finite-size scaling equation ( 126]) . The numerical 
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accuracy one reaches depends highly on the system size reached in the calculation. 

For the 3-12 lattice and the (lxl) : (1 x 1) (the kagome) and (2x2) : (2 x 2) kagome-type 
subnet lattices, the largest dimension of the arrays used to store the values and the positions 
of nonzero elements of the sparse matrices is d^Vy The largest system size we reached is 
L = 15 with d { ™ ] = 335897865. For the (3 x 3) : (3 x 3) and (4 x 4) : (4 x 4) kagome-type 
subnet lattices, the sparse matrix T2 in the magnetic sector is further decomposed in two 
rectangular matrices of dimensions <i^+i x d^7 2 an d ^£+2 x ^l+h an d the largest system size 
we reached is L = 14. The computer memory requirement for the calculations of the largest 
system is about 65 gigabytes, which is quite large, but the CPU time consumed is rather 
modest. It is just a few hours for a typical calculation of the magnetic scaled gap. 

The magnetic scaling dimension is estimated by extrapolating the scaled gaps 
Xh(v c (L), L) at the solution of fl25|) for a sequence of increasing system sizes in accordance 
with (j27|) . The free energy density at the estimated critical point is calculated using (|28|) 
and the conformal anomaly c is computed by making use of the finite-size scaling relation 
(1291) . Details of the data fitting procedure are described in [161]. We also checked corrections 
to scaling due to the leading irrelevant field. Take the simple kagome lattice as an example. 
According to the Coulomb gas theory [29], y t2 =4 — 4/(7 with g given in fl3T|) is the second 
leading thermal exponent, which we expect to be a candidate for the leading correction 
exponent y u . For q > 2, we indeed found y u close to y t2 . For q = 0.5, y u is about —2.00(1), 
which is the analytic one. For q = 1.0 and 1.5, we found y u = —1.79(3) and —1.51(2) respec- 
tively, which dominate and overcome the corresponding y t , 2 . For q = 2, the Ising model, we 
obtain y u = —4.00(1). The amplitudes of yt 2 = —4/3 and the analytic y u = —2 corrections 
vanish. This is understandable for lattices with sixfold rotational symmetry. This picture is 
generally true for all (n x n) : (n x n) kagome subnet lattices. 

We summarize in Tables HHVl numerical results of our calculations on the critical point 
v c , conformal anomaly c, and magnetic scaling dimension Xh together with the universality 
predictions of c and Xh- We have also computed v c using the homogeneity assumption and 



list the results. The v c calculation for the kagome lattice extends those of [30] using Monte 
Carlo renormalization group method and finite-site scaling analysis for q = 1,2,3,4. Our 
study extends to non-integer q and offers results with higher accuracy. 

For q = 2, the Ising model, our numerical estimates of the critical threshold agree with 
the exact critical results up to 7 or 8 decimal numbers. This probably indicates the limit 
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TABLE I: Critical properties of the Potts model on the 3-12 lattice. (H = Homogeneity Assump- 
tion, N = Numerical, T = Theoretical universality prediction.) 






V r (H) 


v r (N) 


c (T) 


c (N) 


X h (T) 


X h (N) 


0.5 


2.007916417382387 


2.00788(1) 


-0.445833945 


-0.4458340(1) 


0.082757037 


0.08276(1) 


1.0 


2 852426157798754 


2 8523883(2) 








104166667 


104167m 


1.5 


3 510849695265078 


3 510825f2l 


288024142 


288024fL) 


116778423 


116778(1") 


2.0 


4.073446135573680 


4.0734460(1) 


0.5 


0.500000(1) 


1/8 


0.12500000(1) 


2.5 


4.574927577671523 


4.574952(3) 


0.66584083 


0.66585(1) 


0.130338138 


0.13033(1) 


3.0 


5.033022514872745 


5.033077(3) 


4/5 


0.800(1) 


2/15 


0.13333(1) 


3.5 


5.458234413883058 


5.458313(2) 


0.910294591 


0.91(1) 


0.133771753 


0.1339(3) 


4.0 


5.857394827983647 


5.857497(3) 


1 


0.999(1) 


1/8 


0.13(1) 



TABLE II: Critical properties of the Potts model on the kagome (1 x 1) lattice. (H = Homogeneity 
Assumption, N = Numerical, T = Theoretical universality prediction.) 



Q 


v c (H) 


v c (N) 


c (T) 


c (N) 


X h (T) 


X h (N) 


0.5 


0.787417375457453 


0.787320(1) 


-0.445833945 


-0.445834(1) 


0.082757037 


0.082757(1) 


1.0 


1.102738621067509 


1.10262924(2) 








0.104166667 


0.104167(1) 


1.5 


1.342082948593078 


1.3420126(2) 


0.288024142 


0.2880243(3) 


0.116778423 


0.116780(3) 


2.0 


1.542459756837412 


1.5424598(1) 


0.5 


0.500000(1) 


1/8 


0.12500000(1) 


2.5 


1.718102046569530 


1.718191(3) 


0.66584083 


0.66584(1) 


0.130338138 


0.1304(1) 


3.0 


1.876269208345760 


1.876458(3) 


4/5 


0.8000(1) 


2/15 


0.1333(1) 


3.5 


2.021253955272383 


2.02154(2) 


0.910294591 


0.910(1) 


0.133771753 


0.134(1) 


4.0 


2.155842236513638 


2.15620(5) 


1 


1.00(1) 


1/8 


0.13(1) 



of the numerical accuracy of the finite- size analysis we can reach at present. For q ^ 2, 
the critical points obtained from fl6]) under the homogeneity assumption coincide with our 
numerical estimations to 5 or so decimal places but lie outside error bars. This indicates that 
the homogeneity assumption, while highly accurate, is an excellent approximation yielding 
numerical values with an error within one part in 10 5 . Our computed values of the conformal 
anomaly c and the magnetic scaling dimension coincide with the theoretical universality 
predictions within error bars. 
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TABLE III: Critical properties of the Potts model on the (2 x 2) : (2 x 2) kagome lattice. (H 
Homogeneity Assumption, N = Numerical, T = Theoretical universality prediction.) 






V r (H) 


Vr (N) 


c (T) 


c (N) 


X h (T) 


X h (N) 


0.5 


1.115482279992555 


1 1154309(3) 


-0 445833945 

\J • J- -L v7 v7 f . J '.J \ ' -L '7 


-0 4458340(2) 


082757037 

V7 • V7 ^7 1 '.7 1 V7 1 


082757(1) 


1.0 


1 505450910604828 

X ■ v7 V7 "7 i v7 V7 ■./ lUUU^UiU 


1 5053987(1) 

-L • '.7 V7 '.7 ? 7 ' 7 V.7 1 1 i J 


o 


o 


104166667 

V7 • -L V / X 1 V7V7 V7Y7 1 


104167(1) 


1.5 


1 790803965420646 


1 7907720(2) 


288024142 


288024(1) 


116778423 


11678(4) 


2.0 


2.024382957091806 


2.02438295(3) 


0.5 


0.500000(1) 


1/8 


0.1250000(1) 


2.5 


2.225885325024986 


2.2259229(2) 


0.66584083 


0.66584(1) 


0.130338138 


0.13034(1) 


3.0 


2.405138877193783 


2.4052181(3) 


4/5 


0.8001(1) 


2/15 


0.1333(1) 


3.5 


2.567855953492942 


2.567981(2) 


0.910294591 


0.910(1) 


0.133771753 


0.1339(3) 


4.0 


2.717691692682905 


2.717856(2) 


1 


0.99(1) 


1/8 


0.13(1) 



TABLE IV: Critical properties of the Potts model on the (3 x 3) : (3 x 3) kagome lattice. (H 
Homogeneity Assumption, N = Numerical, T = Theoretical universality prediction.) 



q 


v c (H) 


v c (N) 


c (T) 


c (N) 


Xh (T) 


X h (N) 


0.5 


1.236699591471530 


1.2366855(3) 


-0.445833945 


-0.4458340(2) 


0.082757037 


0.0827569(2) 


1.0 


1.626971272019731 


1.6269594(2) 








0.104166667 


0.104167(1) 


1.5 


1.906766682469675 


1.906760(1) 


0.288024142 


0.288024(1) 


0.116778423 


0.116779(2) 


2.0 


2.133002727374153 


2.13300273(1) 


0.5 


0.500000(1) 


1/8 


0.1250000(1) 


2.5 


2.326449318777172 


2.32645568(5) 


0.66584083 


0.66585(1) 


0.130338138 


0.130338(1) 


3.0 


2.497336478778200 


2.4973486(2) 


4/5 


0.800(1) 


2/15 


0.1333(1) 


3.5 


2.651556985414795 


2.651575(3) 


0.910294591 


0.91(1) 


0.133771753 


0.1338(1) 


4.0 


2.79285603450327 


2.79288(2) 


1 


0.999(1) 


1/8 


0.13(1) 



Our numerical results for bond percolation are summarized in Table IVll for the 3-12 lattice 
and the (n x n) : (n x n) subnet kagome-type lattices. For the kagome lattice, we found 
p c = 0.524404978(5), which coincides with the best estimation {24]. For the 3-12 lattice our 
numerical result of p c = 0.74042077(2) is in agreement with other findings [ll, Q, 34 1 to 6 
decimal places. For kagome-type subnet lattices, our numerical analysis determines p c with 
an accuracy up to 7 or 8 decimal places. 

In Table IVll we also give thresholds computed using the homogeneity assumption ([6]). The 
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polynomial equations determining the bond percolation thresholds p c under the homogeneity 
assumption ([6]) for (n x n) : (n x n) subnet lattices in Table I VI I are as follows: 

1 - 3p 2 - 6p 3 + I2p 4 - 6p 5 + p 6 = 0, (n = 1), (32) 
1 - 3p 4 - 18p 5 - 39p 6 + 30p 7 + 273/ + 264p 9 - 1785j9 10 - 126p n + 8232p 12 
-162326p 13 + 16359p 14 - 9948p 15 + 3708p 16 - 786p 17 + 73p 18 = 0, (n = 2), (33) 

1 - 3p 6 - 36p 7 - 186/ - 372p 9 + 447p 10 + 3558j9 n + 471 lp 12 - 5274p 13 - 30771p 14 

- 110816J9 15 + 69828p 16 + 1309302j9 17 - 242760p 18 - 10117626p 19 + 9190737j9 20 
+ 53446600p 21 - 137597577p 22 - 15101358p 23 + 714425889p 24 - 1897059306p 25 
+ 2985201585p 26 - 3337272356p 27 + 2817156177p 28 - 1840940730p 29 

+ 938230487p 30 - 371179194p 31 + 112125462p 32 - 25052124p 33 + 3909120p 34 

- 380880J9 35 + 17464p 36 = 0, (n = 3) (34) 



TABLE V: Critical properties of the Potts model on the (4 x 4) : (4 x 4) kagome lattice. (H = 
Homogeneity Assumption, N = Numerical, T = Theoretical universality prediction.) 



q 


v c (H) 


v c (N) 


c (T) 


c (N) 


X h (T) 


X h (N) 


0.5 


1.287715536704650 


1.2877116(2) 


-0.445833945 


-0.4458340(1) 


0.082757037 


0.0827569(1) 


1.0 


1.669262339202358 


1.6692593(3) 








0.104166667 


0.10417(1) 


1.5 


1.941284616762751 


1.9412832(5) 


0.288024142 


0.288024(1) 


0.116778423 


0.11678(1) 


2.0 


2.160721132019555 


2.160721132(1) 


0.5 


0.500000(1) 


1/8 


0.1250000(1) 


2.5 


2.348099505779181 


2.3481001(2) 


0.66584083 


0.66585(1) 


0.130338138 


0.13034(1) 


3.0 


2.513467694176093 


2.5134684(2) 


4/5 


0.800(1) 


2/15 


0.1333(1) 


3.5 


2.662592230189568 


2.662594(3) 


0.910294591 


0.911(1) 


0.133771753 


0.134(1) 


4.0 


2.799129506399588 


2.799132(5) 


1 


0.999(1) 


1/8 


0.13(1) 
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1 - 3p 8 - QOp 9 - 528p w - 2406p n - 4518p 12 + 8388p 13 + 64323p 14 + 108744p 15 - 149520p 16 

- 892404p 17 - 664532p 18 + 2272086p 19 - 2348817p 20 - 12425874p 21 + 123063933p 22 

+ 344663478p 23 - 1382031989p 24 - 5244471786p 25 + 12598666671p 26 + 50539880448p 27 

- 112896871341p 28 - 350902330710p 29 + 955575283123p 30 + 1782743557128p 31 

- 7239409905561p 32 - 576723 1526534p 33 + 52365034246041p 34 - 23401813013430p 35 

- 353073527306441p 36 + 1041090144149322p 37 - 623756767891383j9 38 - 43674772479 15326p 39 
+ 181176071 72859264p 40 - 42034422047996604p 41 + 71675099615055545p 42 

- 97479081216503664p 43 + 109775262989475858p 44 - 104474772230850020p 45 
+ 85036825023972936p 46 - 59604466077733650p 47 + 36101308809040333p 48 

- 18909591474961260p 49 + 8552494666923729p 50 - 3327421649714158p 51 + 1106659102637175p 52 

- 311767535257674p 53 + 73442507365712p 54 - 14206464131418p 55 + 2198697552561p 56 

- 261883431344p 57 + 22544948382p 58 - 1248899580p 59 + 33437353/° = 0, (n = 4). (35) 



The threshold ( 132]) for n = 1 has previously been given in [30, |37| and in 



33 



Thresholds 



331) - ( 1351) for n = 2,3,4 are new. The polynomial equation determining p c under the 



homogeneity assumption in Table |Vl]for the 3-12 lattice has been given in I and 
is 1 — p + p 2 + p 3 — 7p 4 + 4p 5 = 0. 



and 



TABLE VI: The thresholds p c of bond percolation on the 3-12 and various kagome-type subnet 
lattices. (H = Homogeneity Assumption, N = Numerical.) 



Subnet 


Pc (H) 


Pc (N) 


Other sources 


(1 x 1) : (1 x 1) 


0.524429717521274 


0.524404978(5) 


0.52440499(2) [24] 
0.52440503(5) [24] 
0.5244053(3) [32] 


(2 x 2) : (2 x 2) 


0.600870248238631 


0.60086193(3) 


0.6008624(10) [11] 


(3 x 3) : (3 x 3) 


0.619333484666866 


0.61933176(5) 


0.6193296(10) [11] 


(4 x 4) : (4 x 4) 


0.625364661497144 


0.62536424(7) 


0.625365(3) [11] 


3-12 lattice 


0.740423317919897 


0.74042077(2) 


0.74042118 [33] 
0.74042081 [11] 
0.74042195(80) [34] 
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B. Site percolation on the 3-12 lattice 



The exact critical threshold for site percolation on the 3-12 lattice is known to be s c = 
y/l — 2 sin(7r/18). It was first given in 39] and is shown in I to be the same as that of 
the (2 x 2) : (2 x 2) Potts subnet lattice with pure 3-site interactions. To calibrate our 
numerical approach, we have also computed s c using the transfer matrix approach. Our 
numerical determination of critical properties of site percolation is summarized in the last 
row in Table I VIII The comparison of numerical estimates of thresholds with exact results 
shows agreements up to 7 decimal places, indicating our numerical estimates to be accurate 
to the same degree of accuracy. Our numerical determination of the conformal anomaly and 
magnetic scaling dimension of site percolation indicates that these models all belong to the 
two-dimensional q — 1 Potts model universality class. Again, the hypothesis of universality 
is verified. 



TABLE VII: Critical properties of site percolation on (1 x 1) : (n x n) kagome-type subnet lattices 
and the 3-12 lattice. (E = Exact result, N = Numerical, T = Theoretical universality prediction.) 



Subnet 


s c (E) 


s c (N) 


Other sources 


c(T) 


c(N) 


Xh (T) 


X h (N) 


(lxl): (lxl) 


0.652703644666139 


0.6527035(2) 


1 - 2sin(vr/18) [35J 








0.1041667 


0.1042(1) 


(1 x 1) : (2 x 2) 


0.707106781186548 


0.7071068(2) 


1/V2 [36J 








0.1041667 


0.10416(1) 


(1 x 1) : (3 x 3) 


0.728355596425196 


0.7283555(1) 










0.1041667 


0.10417(1) 


(1 x 1) : (4 x 4) 


0.738348473943256 


0.7383483(5) 










0.1041667 


0.10417(2) 


(1 x 1) : (5 x 5) 


0.743548682503071 


0.7435486(3) 










0.1041667 


0.1042(1) 


(1 x 1) : (6 x 6) 


0.746418147634282 


0.7464180(3) 










0.1041667 


0.10417(1) 


















3-12 lattice 


a/1 -2sin(7r/18) 


0.8079008(3) 


Vl-2sm(7r/18) [39] 








0.1041667 


0.10416(1) 




=0.807900764120 




0.807904(4) [N] [39] 











C. Exact thresholds for site percolation on the (1 x 1) : (n x n) kagome-type lattice 

It was shown in I, and in Sec. [TTl that the rigorous critical frontier (Tjf]) yields the exact 
thresholds of site percolation on (1 x 1) : (nxn) kagome-type subnet lattices. The polynomial 
equations determining the threshold s c are generated by substituting expressions of A and C 
in (A. 5) - (A. 11) into (jlj) and setting q — 1. This yields the site percolation thresholds shown 
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in Table IVlTI Explicitly, the thresholds for site percolation on (1 x 1) : (n x n) kagome-type 
subnet lattices, 1 < n < 6, in Table lyTTI are as follows: 

1 - 3s 2 + s 3 = 0, (n = 1) (36) 
1 - 3s 3 - 3s 4 + 6s 5 - 2s 6 = 0, (n = 2) (37) 
l-3s 4 -9s 5 + 9s 6 + lls 7 -12s 8 + s 9 + s 10 = 0, (n = 3) (38) 
1 - 3s 5 - 18s 6 + 12s 7 + 12s 8 + 41s 9 - 66s 10 + 9s 11 - 9s 12 + 48s 13 - 36s 14 + 8s 15 = 0, 

(n = 4) (39) 
1 - 3s 6 - 30s 7 + 60s 9 - 30s 10 + 216s 11 - 329s 12 - 48s 13 + 3s 14 + 396s 15 

+180s 16 - 1113s 17 + 1038s 18 - 393s 19 + 48s 20 + 3s 21 = 0, (n = 5) (40) 
1 - 3s 7 - 45s 8 - 45s 9 + 165s 10 - 75s 11 + 165s 12 + 510s 13 - 1056s 14 - 959s 15 + 367s 16 
+2349s 17 + 3433s 18 - 6589s 19 - 9069s 20 + 22070s 21 - 11495s 22 - 3597s 23 + 4455s 24 
+702s 25 - 1971s 26 + 792s 27 - 106s 28 = 0, (n = 6). (41) 



The threshold s c = 1 — 2sin(7r/18) for the n = 1 kagome lattice was first given in 35] . The 



;hreshold s c — 1/ \/2 for n = 2 has also been obtained by a "cell-to-cell" transformation in 



361 ] . Here, the thresholds for 3 < n < 6 are new. 

We also computed s c and other critical properties numerically. The results are summa- 
rized in Table IVHI Again, our numerical estimates of s c agree with the exact results up to 
7 decimal places, and these models all belong to the two-dimensional q = 1 Potts model 
universality class. 

Finally, we comment on some numerical specifics. Since the number of non-crossing 
partitions for site percolation is much larger than that of the Potts model in both the 
magnetic and non-magnetic sector for a given circumference L, the maximum system size 
L = 12 that we reached is smaller. The largest dimension of arrays used to save values and 
positions of nonzero elements of the sparse matrices is df^ 3 = 125481607, which requires 
about 43 gigabytes computer memory in the calculations. Corrections to scaling due to the 
leading irrelevant field is about —1.8(1) for all lattices. 



21 



VI. SUMMARY 

We have studied critical properties of the g-state Potts model and bond and site perco- 
lation on two general classes of lattices, the triangular-type and kagome-type lattices. For 
the triangular- type lattices of Fig. 1(a), the exact critical frontier is known and this led 
to a determination of the exact critical thresholds of site percolation on (1 x 1) : (n x n) 
kagome-type subnet lattices. Results for 1 < n < 6 are given. For the kagome-type lattices 
of Fig. 1(b), no exact results are known except for q = 2. We carried out finite-size analysis 
to numerically determine critical properties for various lattice models including the 3-12 lat- 
tice and (n x n) : (n x n) kagome-type subnet lattices. Our numerical results on conformal 
anomaly and magnetic correlation length verify that the principle of universality holds. 

We have also computed the critical thresholds for the Potts and bond percolation on the 3- 
12 lattice and (nxn) : (nxn) kagome-type subnet lattices using the homogeneity assumption 
OH]). To assess the accuracy of our numerical analysis as well as that of the homogeneity 
assumption, we have applied our numerical procedure to study critical properties of models 
for which exact results are known. The comparison of numerical and known results shows 
that the numerical procedure is accurate to 7 or 8 significant digits in determining critical 
thresholds. Assuming the same degree of accuracy for all lattices, this in turn infers that 
the homogeneity assumption determines critical threshold with an accuracy up to 5 decimal 
places or higher. 

Finally, our analysis of critical properties is based on the use of lattice-dependent con- 
stants A, B, C for the hatched triangles shown in Fig. 1. We have developed an algorithm 
of evaluating expressions of A,B,C using computers for hatched triangles in the form of a 
stack-of-triangle structure. 
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Appendix: Expressions of A, B, C 

In this Appendix we list constants A, B, C computed using the computer algorithm as 
described in Sec. III. The condition ([3]) holds for for ferromagnetic Potts models with 
v,m,q > 0. This confirms that qA = C is the exact critical frontier in the ferromagnetic 
regime. 

1. Potts model on n x n subnets with pure 2-site coupling K 

v = e K - 1 
Subnet lxl: 

A=l, B = v, C = v 3 + 3v 2 . (A.l) 

Subnet 2x2: 

A = q 3 + 9q 2 v + 33qv 2 + (50 + Aq)v 3 + 2lv 4 + 3v 5 
B = q 2 v 2 + 10qv 3 + (30 + 2q)v 4 + 22v 5 + 7v 6 + v 7 
C = 9qv 4 + (54 + 3q)v 5 + 63v 6 + 33v 7 + 9v 8 + v 9 . (A.2) 

Subnet 3x3: 

A = 29q 3 v 6 + (459g 2 + 9q 3 )v 7 + (2592g + 423q 2 )v 8 + (5292 + 4185g + 171q 2 )v 9 
+ (12825 + 3258g + 36q 2 )v 10 + (15534 + 1539g + 3q 2 )v u + (12184 + 454g)t» 12 
+ (6732 + 78q)v 13 + (2688 + 6q)v u + 768v 15 + I50v w + I8v 17 + v 18 , 

B = q b v 3 + 21gV + (199g 3 + 3q 4 )v 5 + (1040g 2 + 85q 3 )v 6 + (2979g + 844g 2 
+ 17q 3 )v 7 + (3780 + 3834g + 340g 2 + 2q 3 )v 8 + (7182 + 2429g + 86q 2 )v 9 
+ (6858 + 950g + 13q 2 )v 10 + (4250 + 233g + q 2 )v u + (1846 + 33q)v 12 
+ (570 + 2q)v 13 + 121v 14 + 16w 15 + v w , 

C = 29q 3 v 6 + (459g 2 + 9q 3 )v 7 + (2592g + 423g> 8 + (5292 + 4185g + 171g 2 )t» 9 
+ (12825 + 3258g + 36g 2 )i; 10 + (15534 + 1539g + 3q 2 )v 11 + (12184 + AhAq)v 12 
+ (6732 + 78q)v 13 + (2688 + 6q)v 14 + 768v 15 + I50v 16 + I8v 17 + v 18 . (A.3) 
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Subnet 4x4: 



A = q 12 + 30q u v + 435g 1( V + (4044g 9 + lQq 10 )v 3 + (26952g 8 + 450g> 4 + (136323g 7 
+ 6057g 8 + 18g> 5 + (539687g 6 + 51643g 7 + 606g 8 )w 6 + (1696074g 5 + 310764g 6 
+ 9336g 7 + 27q 8 )v 7 + (4229307g 4 + 1388580g 5 + 88062g 6 + 9A8q 7 )v 8 + (8218405g 3 
+ 4700302g 4 + 568443g 5 + 15274g 6 + A3q 7 )v 9 + (11888619g 2 + 11989938g 3 
+ 2635653g 4 + 148602g 5 + 1608g 6 )w 10 + (11554638g + 22208904g 2 + 8874255g 3 
+ 966189g 4 + 26310g 5 + 96q 6 )v n + (5728860 + 27132098g + 21124723g 2 + 4354609g 3 
+ 251967g 4 + 3222g 5 + 3q 6 )v 12 + (16691130 + 32535501g + 13440978g 2 + 1559913g 3 
+ 48486g 4 + 2Q7q 5 )v 13 + (24925347 + 26263347g + 6337347g 2 + 426225g 3 + 6834g 4 
+ I2q 5 )v 14 + (25218686 + 15860794g + 2323668g 2 + 89613g 3 + 644g> 15 
+ (19264962 + 7537923g + 674634g 2 + 14088g 3 + 30q 4 )v 16 + (11718378 + 2882733g 
+ 154248g 2 + 1539g 3 )v 17 + (5825765 + 890008g + 27036g 2 + 102g 3 )v 18 + (2389554 
+ 219258g + 3447g 2 + 3q 3 )v 19 + (806778 + 42009g + 288g 2 )v 20 + (221570 + 5996g 
+ 12q 2 )v 21 + (48435 + 594q)v 22 + (8136 + 36q)v 23 + (990 + q)v 24 + 78v 25 + 3v 26 , 

B = q 9 v 4 + 36gV + (609g 7 + 4g> 6 + (6340g 6 + 193g> 7 + (44883g 5 + 3730g 6 
+ 30q 7 )v 8 + (225145g 4 + 40957g 5 + 1156g 6 + 3q 7 )v 9 + (804174g 3 + 290085g 4 
+ 19521g 5 + 240g> 10 + (1980543g 2 + 1384197g 3 + 192846g 4 + 6208g 5 + 30g> n 
+ (3064302g + 4394962g 2 + 1220478g 3 + 85315g 4 + 1355g 5 + 2q 6 )v 12 + (2280420 
+ 8560443g + 4979114g 2 + 717391g 3 + 26906g 4 + 201q 5 )v 13 + (7901226 + 12201635g 
+ 3771317g 2 + 307105g 3 + 6145g 4 + 19q 5 )v 14 + (14006718 + 11666444g + 2106366g 2 
+ 99167g 3 + 989g 4 + q 5 )v 15 + (16765996 + 8289847g + 906964g 2 + 24293g 3 
+ 102g> 16 + (15077600 + 4592110g + 306788g 2 + 4431g 3 + 5q 4 )v 17 + (10735261 
+ 2026213g + 81784g 2 + 574g 3 )w 18 + (6217796 + 715829g + 17041g 2 + A8q 3 )v 19 
+ (2966339 + 200961g + 2720g 2 + 2q 3 )v 20 + (1168690 + 43935g + 319g 2 )w 21 
+ (378274 + 7210g + 25q 2 )v 22 + (99306 + 833g + q 2 )v 23 + (20684 + 60q)v 24 
+ (3298 + 2q)v 25 + 379v 26 + 28v 27 + v 28 , 
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C = 99gV + (2871g 5 + 29g> 9 + (36285g 4 + 2250q 5 )v 10 + (256785g 3 + 47886g 4 

+ 765<T> n + (1078677g 2 + 489164g 3 + 30525g 4 + 135q 5 )v 12 + (2567538g 

+ 2747088g 2 + 463617g 3 + 11898g 4 + 9q 5 )v 13 + (2723220 + 8371674g 

+ 3520566g 2 + 282843g 3 + 2919g> 14 + (11070162 + 13828185g + 2965938g 2 

+ 120355g 3 + 420g 4 )w 15 + (22921893 + 15205050g + 1805958g 2 + 36504g 3 + 27q 4 )v w 

+ (31898508 + 12327042g + 826827g 2 + 7773g 3 )t» 17 + (33199952 + 7729828g + 287940g 2 

+ 1098g 3 y 8 + (27253662 + 3833517g + 75627g 2 + 90q 3 )v 19 + (18157536 + 1513779g 

+ 14556g 2 + 3q 3 )v 20 + (9965342 + 473460g + 1938g> 21 + (4531923 + 115287g 

+ 159g 2 )w 22 + (1706052 + 21150g + 6q 2 )v 23 + (527795 + 2757g)v 24 + (132300 + 228q)v 25 

+ (26256 + 9q)v 26 + 3976w 27 + 432w 28 + 30v 29 + v 30 (A.4) 



2. Potts model on n x n subnets with pure 3-site coupling M 



rn 



Subnet lxl: 



A=l, 



B = 0, C = m. 



(A.5) 



Subnet 2x2: 



A = q 3 + 3qm, B = m 2 , C = m 3 . 



(A-6) 



Subnet 3x3: 



A = q 7 + Qq 5 m + 15g 3 m 2 + (Uq + 3q 2 )m 3 + 3m 4 , 
B = q 2 m 3 + (2 + 2g)m 4 + m 5 , 
C = 3m 5 + m 6 . 



(A.7) 
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Subnet 4x4: 

A = q 12 + 10g 10 m + 45g 8 m 2 + (114g 6 + 6g 7 )m 3 + (165g 4 + 42g 5 )m 4 

+ (117g 2 + 99g 3 + 9g 4 )m 5 + (20 + 73q + 33g 2 + 3g 3 )m 6 

+ (15 + 13g)m 7 + 3m 8 , 
B = q 5 m 4 + (6g 3 + 3g 4 )m 5 + (9g + 15g 2 + 3g 3 )m 6 + (12 + 12g + 3g 2 )m 7 

+ (6 + 2g)m 8 + m 9 

C = (2 + 9g)m 7 + (15 + 3g)m 8 + 7m 9 + m 10 . (A.8) 
Subnet 5x5: 

A = q 18 + 15g 16 m + 105g 14 m 2 + (445g 12 + 10g 13 )m 3 + (1245g 10 + 120g n )m 4 + (2358g 8 
+ 624g 9 + 18g 10 )m 5 + (2967g 6 + 1795g 7 + 189g 8 + 6g 9 )m 6 + (2298g 4 + 2976g 5 + 792g 6 
+ 84g 7 )m 7 + (888g 2 + 2613g 3 + 1608g 4 + 399g 5 + 27g 6 )m 8 + (86 + 864g + 1416g 2 
+ 808g 3 + 182g 4 + 9g 5 )m 9 + (249 + 579g + 345g 2 + 81g 3 + 3g 4 )m 10 + (126 + 129g 
+ 36g 2 )m n + (30 + 13g)m 12 + 3m 13 , 

B = q 9 m 5 + (12g 7 + 4g 8 )m 6 + (47g 5 + 42g 6 + 6g 7 )m 7 + (80g 3 + 156g 4 + 57g 5 + 7g 6 )m 8 
+ (42q + 221g 2 + 204g 3 + 59g 4 + 7g 5 )m 9 + (67 + 225g + 169g 2 + 58g 3 + 5g 4 )m 10 
+ (105 + 121g + 37g 2 + 5g 3 )m n + (49 + 25g + 3g 2 )m 12 + (11 + 2g)m 13 + m 14 , 

C = (12q 2 + 29g 3 )m 9 + (6 + 63g + 96g 2 + 9g 3 )m 10 + (102 + 132g + 36g 2 )m n + (127 + 51g 
+ 3g 2 )m 12 + (57 + 6q)m 13 + 12m 14 + m 15 . (A.9) 



26 



Subnet 6x6: 

A = q 25 + 21g 23 m + 210g 21 m 2 + (1315g 19 + 15g 20 )m 3 + (5715g 17 + 270g 18 )m 4 

+ (18084g 15 + 2235g 16 + 30g 17 )m 5 + (42560g 13 + 11166g 14 + 525g 15 + 10g 16 )m 6 

+ (74769g n + 37113g 12 + 4113g 13 + 210g 14 )m 7 + (96819g 9 + 85161g 10 + 18795g n 

+ 1926g 12 + 54g 13 )m 8 + (89350g 7 + 135082g 8 + 54561g 9 + 10086g 10 + 798g n 

+ 18g 12 )m 9 + (54795g 5 + 143112g 6 + 102321g 7 + 32634g 8 + 5172g 9 + 351g 10 

+ 6g n )m 10 + (19257g 3 + 92445g 4 + 118923g 5 + 65685g 6 + 18570g 7 + 2622g 8 

+ 174g 9 )m n + (2642g + 29213g 2 + 74961g 3 + 76633g 4 + 38522g 5 + 9990g 6 + 1452g 7 

+ 66g 8 )m 12 + (2106 + 17556g + 41547g 2 + 42069g 3 + 20271g 4 + 5523g 5 + 735g 6 

+ 27g 7 )m 13 + (4536 + 16626g + 18789g 2 + 10158g 3 + 2817g 4 + 381g 5 + 9g 6 )m 14 

+ (3459 + 6819g + 4372g 2 + 1398g 3 + 179g 4 + 3g 5 )m 15 + (1398 + 1593g + 594g 2 

+ 87g 3 )m 16 + (342 + 216g + 36g 2 )m 17 + (48 + 13g)m 18 + 3m 19 , 

B = g 14 m 6 + (20g 12 + 5g 13 )m 7 + (145g 10 + 90g n + 10g 12 )m 8 + (540g 8 

+ 623g 9 + 168g 10 + 14g n )m 9 + (1142g 6 + 2263g 7 + 1132g 8 + 221g 9 + 17g 10 )m 10 
+ (1304g 4 + 4549g 5 + 4089g 6 + 1413g 7 + 258g 8 + 16g 9 )m n + (659g 2 + 4573g 3 
+ 7808g 4 + 4873g 5 + 1569g 6 + 238g 7 + 16g 8 )m 12 + (68 + 1652g + 6525g 2 + 8325g 3 
+ 4896g 4 + 1432g 5 + 214g 6 + 14g 7 )m 13 + (1183 + 4981g + 6786g 2 + 3999g 3 + 1184g 4 
+ 182g 5 + llg 6 )m 14 + (1950 + 4062g + 2733g 2 + 885g 3 + Ulq 4 + 9g 5 )m 15 + (1372 
+ 1566g + 571g 2 + 102g 3 + 7g 4 )m 16 + (538 + 338g + 64g 2 + 5g 3 )m 17 + (127 + 40g 
+ 3g 2 )m 18 + (17 + 2g)m 19 + m 20 , 

C = (54g 5 + 99g 6 )m n + (146g 3 + 585g 4 + 483g 5 + 29g 6 )m 12 + (90g + 888g 2 
+ 1971g 3 + 1170g 4 + 159g 5 )m 13 + (228 + 1992g + 3303g 2 + 2073g 3 + 381g 4 
+ 9g 5 )m 14 + (1677 + 3726g + 2598g 2 + 660g 3 + 33g 4 )m 15 + (2421 + 2511? 
+ 810g 2 + 81g 3 )m 16 + (1626 + 813g + 114g 2 + 3g 3 )m 17 + (620 + 133g + 6g 2 )m 18 
+ (141 + 9g)m 19 + 18m 20 + m 21 . (A.10) 
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Subnet 7x7: 



A = g 33 + 28g 31 m + 378g 29 m 2 + (3255g 27 + 21g 28 )m 3 + (19950g 25 

+ 525g 26 )m 4 + (92025g 23 + 6210g 24 + 45g 25 )m 5 + (329615g 21 + 45955g 22 

+ 1155g 23 + 15g 24 )m 6 + (932792g 19 + 237060g 20 + 13755g 21 + 430g 22 )m 7 

+ (2102508g 17 + 899163g 18 + 100401g 19 + 5835g 20 + 90g 21 )m 8 + (3777012g 15 

+ 2577591g 16 + 499731g 17 + 48899g 18 + 2062g 19 + 30g 20 )m 9 + (5372928g 13 

+ 5650464g 14 + 1784592g 15 + 278639g 16 + 22401g 17 + 846g 18 + 10g 19 )m 10 

+ (5963664g n + 9467118g 12 + 4680099g 13 + 1127826g 14 + 150330g 15 + 10692g 16 

+ 396g 17 )m n + (5030270g 9 + 11966359g 10 + 9054966g n + 3301756g 12 

+ 683417g 13 + 80559g 14 + 5754g 15 + 132g 16 )m 12 + (3081149g 7 + 11079301g 8 

+ 12758268g 9 + 6987391g 10 + 2179863g n + 399888g 12 + 46751g 13 + 2766g 14 

+ 54g 15 )m 13 + (1263348g 5 + 7101903g 6 + 12629796g 7 + 10471011g 8 

+ 4893192g 9 + 1358682g 10 + 242670g n + 26103g 12 + 1422g 13 + 18g 14 )m 14 

+ (294802g 3 + 2824768g 4 + 8154549g 5 + 10567810g 6 + 7525012g 7 + 3161463g 8 

+ 844184g 9 + 144491g 10 + 15006g n + 690g 12 + 6g 13 )m 15 + (25948g + 547184g 2 

+ 2935194g 3 + 6464096g 4 + 7397804g 5 + 4847082g 6 + 1962007g 7 + 511413g 8 

+ 86643g 9 + 8408g 10 + 348g n )m 16 + (23688 + 392541g + 1878534g 2 + 3986040g 3 

+ 4432011g 4 + 2878479g 5 + 1160562g 6 + 304893g 7 + 50895g 8 + 4761g 9 

+ 153g 10 )m 17 + (112270 + 781871g + 1899615g 2 + 2296835g 3 + 1561921g 4 

+ 654696g 5 + 175858g 6 + 29445g 7 + 2616g 8 + 66g 9 )m 18 + (162018 + 662861g 

+ 987349g 2 + 764529g 3 + 345520g 4 + 97287g 5 + 16524g 6 + 1404g 7 + 27g 8 )m 19 

+ (123912 + 327444g + 320334g 2 + 166830g 3 + 5081 lg 4 + 8838g 5 + 747g 6 

+ 9g 7 )m 20 + (60508 + 106552g + 69966g 2 + 24288g 3 + 4547g 4 + 372g 5 

+ 3g 6 )m 21 + (20268 + 23740g + 10182g 2 + 2142g 3 + 188g 4 )m 22 + (47 1 

+ 3525g + 894g 2 + 87g 3 )m 23 + (732 + 316g + 36g 2 )m 24 + (69 + 13g)m 25 + 3m 26 , 
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B = q 20 m 7 + (30g 18 + 6g 19 )m 8 + (345g 16 + 165g 17 + 15g 18 )m 9 

+ (2175g 14 + 1820g 15 + 390g 16 + 25g 17 )m 10 + (8587g 12 + 11240g 13 + 4199g 14 

+ 624g 15 + 35g 16 )m n + (22421g 10 + 43891g n + 25772g 12 + 6546g 13 + 844g 14 

+ 40g 15 )m 12 + (38966g 8 + 113011g 9 + 100896g 10 + 39578g n + 8607g 12 + 951g 13 

+ 43g 14 )m 13 + (43431g 6 + 190110g 7 + 258498g 8 + 153642g 9 + 50665g 10 + 9597g n 

+ 981g 12 + 44g 13 )m 14 + (28116g 4 + 196731g 5 + 420318g 6 + 385688g 7 + 190950g 8 

+ 55686g 9 + 9626g 10 + 970g n + 40g 12 )m 15 + (8320g 2 + 108004g 3 + 395182g 4 

+ 590541g 5 + 455390g 6 + 203444g 7 + 54523g 8 + 9175g 9 + 874g 10 + 37g n )m 16 

+ (488 + 20976g + 170827g 2 + 477931g 3 + 626376g 4 + 452018g 5 + 190124g 6 + 50024g 7 

+ 8064g 8 + 775g 9 + 32g 10 )m 17 + (15646 + 137655g + 396335g 2 + 531376g 3 + 381398g 4 

+ 162295g 5 + 42103g 6 + 6904g 7 + 640g 8 + 28g 9 )m 18 + (51536 + 224439g + 352647g 2 

+ 280816g 3 + 123719g 4 + 33731g 5 + 5505g 6 + 532g 7 + 22g 8 )m 19 + (66881 + 176633g 

+ 170578g 2 + 86821g 3 + 24401g 4 + 4339g 5 + 407g 6 + 18g 7 )m 20 + (49197 + 83330g 

+ 51390g 2 + 17058g 3 + 3016g 4 + 326g 5 + 13g 6 )m 21 + (23535 + 25741g + 10127g 2 

+ 2132g 3 + 217g 4 + llg 5 )m 22 + (7727 + 5333g + 1281g 2 + 155g 3 + 7g 4 )m 23 

+ (1751 + 719g + 94g 2 + 5g 3 )m 24 + (264 + 57g + 3g 2 )m 25 + (24 + 2g)m 26 + m 27 , 

C = (222g 9 + 351g 10 )m 13 + (1386g 7 + 3816g 8 + 2250g 9 + 99g 10 )m 14 

+ (3198g 5 + 14696g 6 + 18894g 7 + 7530g 8 + 686g 9 )m 15 + (3072g 3 + 24933g 4 
+ 59031g 5 + 51366g 6 + 18519g 7 + 2208g 8 + 29g 9 )m 16 + (888g + 16710g 2 + 78501g 3 
+ 132330g 4 + 99747g 5 + 35370g 6 + 5022g 7 + 150g 8 )m 17 + (1954 + 32712g + 133203g 2 
+ 202318g 3 + 150051g 4 + 55086g 5 + 8997g 6 + 456g 7 )m 18 + (25370 + 130809g + 228960g 2 
+ 178833g 3 + 72045g 4 + 13566g 5 + 948g 6 + 9g 7 )m 19 + (67692 + 179352g + 174174g 2 
+ 78354g 3 + 17760g 4 + 1581g 5 + 33g 6 )m 20 + (81987 + 127689g + 73398g 2 + 19325g 3 
+ 2343g 4 + 72g 5 )m 21 + (58543 + 54880g + 18612g 2 + 2746g 3 + 150g 4 )m 22 + (27453 
+ 15033g + 2859g 2 + 207g 3 + 3g 4 )m 23 + (8825 + 2607g + 246g 2 + 6g 3 )m 24 + (1949 
+ 264g + 9g 2 )m 25 + (285 + 12g)m 26 + 25m 27 + m 28 . (A.ll) 
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